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Abstract 

The thermodynamics of a free Bose gas with effective temperature scale T and hard-sphere Bose 
gas with the T scale are studied. T arises as the temperature experienced by a single particle in a 
quantum gas with 2-body harmonic oscillator interaction V osc , which at low temperatures is expected 
to simulate, almost correctly, the attractive part of the interatomic potential Vhc between 4 He atoms. 
The repulsive part of Vhc is simulated by a hard-sphere (HS) potential. The thermodynamics of this 
system of HS bosons, with the T temperature scale (HSET), is investigated, first, by the Bogoliubov- 
Huang method and next by a modified version of this method, which takes approximate account of 
those terms of the 2-body repulsion which are linear in the zero-momentum Bose operators ao, ag 
(originally rejected by Bogoliubov). Theoretical heat capacity Cv(T) exhibits good agreement, below 
2.1 K, with the experimental heat capacity graph observed in 4 He at saturated vapour pressure. The 
phase transition to the low-temperature phase , with a Bose-Einstein condensate, occurs in the HSET 
at T\ =2.17 K, and is accompanied, in the modified HSET version, by a singularity of CV(T). Other 
thermal properties of HSET, such as the momentum distribution function, the fraction of atoms in the 
momentum condensate and normal fluid density, agree qualitatively with those of 4 He, but improve 
those of the free Bose gas. 

PACS numbers: 64.70 Dv, 61.12 Ex, 67.80 Gb 



1 Introduction 

The extraordinary low-temperature properties of 4 He have so far evaded full explanation by a first- 
principles microscopic theory, although the fundamental concepts of such theory are well established. 
Neutron scattering experiments from liquid helium [T] have confirmed London's hypothesis [2], which 
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points to Bosc-Einstein condensation (BEC) of 4 Hc atoms below T\ = 2.17 K, as the decisive factor 
behind these properties. Theory shows that a free gas of zero-spin bosons (FBG), with particle mass and 
density equal to those of liquid 4 He, exhibits the BEC transition at To = 3.1486 K. The shift To — T\ > 
is therefore presumably a consequence of interatomic interactions. 

A decisive step in the understanding of BEC of an interacting gas was made by Bogoliubov [3], [4|. 
His ideas were incorporated into the theory of hard-sphere Bose gas (HSBG) by Lee et al. [5] - jS] and 
were refined by Angelescu et al. (§]. Unfortunately, the temperature Tbe at which an interacting Bose 
gas, with 4 He parameters, exhibits BEC, exceeds To (Ref. [7]) or remains equal to To (Ref. [9]) in these 
theories. The HSBG heat capacity in the vicinity of T = is proportionate to T 3 , as in 4 He, [7] and all 
atoms belong to the momentum condensate at T = K. The experimentally estimated fraction of 4 Hc 
atoms in the momentum condensate (FMC) and fraction of 4 He normal fluid (FNF) at T = K are equal 
0.07-0.09 [10] and zero, respectively [IT] . 

A remarkable success was achieved by Ceperley and Pollock [T^] - [H] who computed the canonical 
density matrix for 64 and 125 interacting 4 He atoms, using a Monte Carlo path-integral technique. The 
resulting heat capacity and FNF are in good quantitative agreement with experiment. 

Our objective is to develop a relatively simple theory which would explain, at least some of the unusual 
properties of liquid 4 He, with accuracy comparable to that achieved in Refs. [Jj2] - [2]. This is done 
by resorting to the concepts of HSBG theory (Refs. [5] - [5]), which takes into account the repulsive 
hard-core part of the 2-body interaction Vn e between 4 He atoms. The attractive part of Vn e is simulated 
by a harmonic oscillator 2-body potential V osc , which gives rise to an effective temperature T of the gas. 
T is the effective temperature experienced by a single particle in a quantum Boltzmann gas with 2-body 
oscillator interactions. It is found in Section 2 by examining the structure of the 1 — particle reduced 
density matrix g^ 1 ' of the canonical ./V — particle density matrix of such gas with, additionally, a 
harmonic oscillator attraction Vq between each particle and a fixed centre. For large TV, one finds that, 
for sufficiently weak Vq, g^ w const exp[— /3i?Q ], where = —h 2 A/2m and 

p= (kBT)- 1 = 7- 1 tanh( 7 /3) (1) 

with 

h [Nk s , 

l=- 2 \j-, P = (kBT)-\ (2) 

k denoting the coupling constant of V osc . 

The procedure of reducing q( n ~> amounts (up to a multiplicative constant |A| _JV+1 , where |A| denotes 
the system's volume) to averaging q( n ' over the positions of N — 1 particles of the iV-particle system. 
Thus, from the viewpoint of 1-particle measurements, a quantum Boltzmann gas of particles interacting 
via V^ sc behaves like a gas of free particles with the same mass but at a different temperature T. We 
put forward a conjecture that a Bose gas with the 2-body interaction V osc behaves analogously. This 
assumption is the starting point for further investigation. 
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The next stages of this inquiry consist in deriving the thermodynamics of a free Bose gas with the 
temperature scale T (FBET) (Sec. 3) and hard-sphere Bose gas with this temperature scale (HSET) 
according to the Bogoliubov-Huang approach (Sec. 4) . It is shown that the HSET allows to improve the 
FBET description of low-temperature thermodynamics of 4 He. In particular, the HSET heat capacity 
shows quantitative agreement with the low-temperature branch of experimental 4 He heat capacity below 
1.8 K. Other HSET thermal properties agree qualitatively with those of 4 He, e.g., FMC and FNF, at 
T = OK, equal 0.5224 and 0.1896, respectively. 

An alternative version of the HSET, which takes approximate account of those terms of the 2-body 
repulsion, which are linear in the zero-momentum Bose operators oq, (originally rejected by Bogoliubov 
in Ref. [3]) is studied in Sec. 5. This modified HSET (DHSET) yields singularity of heat capacity at 
T\ and CV(T) agrees with experiment below 2.1 K. The FMC and FNF in this approach equal 0.4266 
and 0.1178, respectively, at T = K. The DHSET thus considerably improves the HSET description of 
helium II thermodynamics. 



2 The effective temperature 



The interatomic potential Vh c has a typical well with a minimum at the distance rg = 3A separating two 
4 He atoms |15j . In the low-temperature regime, a reliable approximation to the partition function 

Z = Tr(cxp[-/3iJ (Ar )]S*W), (3) 

where 

N 

ffW=£F (ri)+ tWIri-rjl) (4) 

i=l l<i<j<N 

and denotes the symmetrizer, can be therefore obtained by applying the Laplace method, i.e., by 
expanding Vn e (r) in the neighbourhood of tq up to second order and integrating the resulting Gaussian 
function. This procedure amounts to replacing Vhc(|i"i — 1*2 1 ) 5 up to an additive constant, by 

^(|r 1 -r 2 |) = ^(|r 1 -r 2 |-r ) 2 , k > 0. (5) 

The approximating partition function 

Z 2 = Tr(exp[-/?fff (6) 

where 

N 

ffW=£tf (r,)+ J2 ^(K-r.-l), (7) 

i=l l<i<j<N 

is intractable, unless we put ro = 0. We are then faced with the problem of solving the thermodynamics 
of TV bosons coupled by harmonic oscillator interactions. The eigenproblem for the Hamiltonian 

H^=H^\ ro=0 (8) 
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of such system, with the Hilbert space L 2 (R 3 ) N , was solved in Ref. [16]. The grand partition 
function Z osc for i?2 ; with ro = 0, expresses in closed form in terms of the fugacity z. However, the 
equation for z takes the same form as for uncoupled oscillators ( Ref. |16j). viz., 

00 

N= Yj (z" 1 cxp[f3fuj(n x + n y + n z )} - ly 1 , (9) 

Tlx ,n y ,n z —0 

and is also intractable, unless approximations are applied (e.g. Ref.|17j). 

The eigenstructure of S^H { 2 N) S^ N \ with r = 0, found in Ref. [H], appears to be inadequate for the 
description of 4 He thermodynamics, also from the viewpoint of existing theories of BEC in homogeneous 
systems, which usually exploit the Hamiltonian +vj i N * > of N bosons, with an interaction Va, enclosed 
in a cube A = L 3 , under periodic boundary conditions. Owing to the different Hilbert spaces of S^ N ^H^ 
and S^H^p, the thermodynamics of the first Hamiltonian does not pass over to that of the later in the 
limit k — > 0, as one would require for a Bose gas with weak V2| ro =o and the usual infinite- volume limit: 
|A| — > 00, N — > 00, iV/| A| = const, cannot be performed for S^ N 'H^[ C . 

A possible approximation to 

e W = ^exp[-/}ffW] (10) 

in S( NS> L 2 (A) N , which fulfils this consistency requirement, can be obtained by exploiting the structure of 
the reduced 1-particle density matrix of the density operator 

g^ = Z^ 1 eM~PH[ N) ], Z x = TreM-PH[ N) ], (11) 

defined in L 2 (R 3 ) N , where 

H[ N) = f>o(r,) + |r?) + \ £ (r.-r,) 2 . (12) 
The reduced density matrices of g^ were studied in Refs. [T5]-[3T]. It was shown in Refs. [TB], [TS] that 

-1 n (N) 



the 1-particle reduced density matrix g^ = Lj^g^ {L X N denoting the partial trace over L 2 (M. 3 ) N 1 ) has 



the form 



gW = Cr i = H « + k -lr 2 , (13) 



where 

.1^ , N -3 



Ci = (2sinh(-Oi)) , O x = tfalfc/m, a\ = ^pmkxjh (14) 



and k, ko, ki, ft, /?i are related by the following formulae 



h = fco4 /2 tanh(^ jv ) + k N kl /2 {N - 1) tanh(foo) 
kH 2 tanh(if2jv) + k^ 2 (N - 1) tanh(±O ) 

n 2 ^fhPi [h^ 

. coth — — \ — 
Vk k N v 2 V to 
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( A^ + (A r_i) ^l + (N-l) 7^=^ , (16) 



with 



Uo = np\ —, n N = hpJ—, k N = k + Nk. (17) 

V m V to 

In the limit /co ~~ ^ 0, fci approaches zero linearly in fco : 

fci = ak + O(fcg) as fc -> 0, (18) 



where 



tanh( 7 /3) + (iV- 1) 7 /3 » ftg x 

tanh( 7 /3) ' 7 ~ 2V m • 1 j 

Using Eqs. (16), (18), one finds 

Eqs. (12), (13), (18), (20) show that, if measurements are restricted to 1-particle observables, a single 
particle of a quantum Boltzmann gas in the state 

gi N) =Z- 1 exp[-/3 J fff >| ro=0 ] (21) 

behaves like one belonging to a free quantum gas at an effective temperature T\ = (A:b/?i(/3)) _1 - Wc 
expect that for a macroscopic sample of a Bose gas with the Hamiltonian (12), the meaning of $\ is 
analogous. In other words, by resorting to the asymptotic form of the 1-particle reduced density matrix 
of a free Bose gas [22], we put forward the following conjecture: for sufficiently small ko > 0, 

Z^LUS W eM-^i N) ])^^(l -*#>)-*, (22) 

where 

Z xs = Tr(SW exphSflf >]), = exp(-/3 1 ij( 1) ). 

Clearly, if T = (fce/3) -1 denotes the real temperature, then T\ > T, if 7 > 0, and 

lim fi = T, (23) 

7— >o 

as required by consistency. 

It is worth noting that the effective temperature of a quantum gas is not an entirely new concept. For 
example Landau and Lifshitz in Ref. |23| (Sec. 33) computed the effective temperature of an interacting 
quantum gas, using a semi-classical approximation. 

Let us now define 

P{fi) = lim ft (/3) = 7 - 1 tanh( 7 /?). (24) 

iV->oo 

This limiting form f3 of /?i provides the best agreement with experimental data on He. Clearly, /? — > /3 
as 7 — > 0. This property of /3, as well as the construction of qY , allow to expect that a sufficiently weak 
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harmonic oscillator interaction, between particles of a quantum gas in equilibrium, can be simulated by 
performing the substitution /3 — > /3. In this approach, the thermodynamics of a gas enclosed in a finite 
volume A, with the 2-body oscillator interaction accounted for by the temperature scale T, passes over 
to the thermodynamics of a free gas in A, with the real scale T, when this interaction is switched off. In 
this manner, the consistency requirement, discussed above, is fulfilled. The consequences of this approach 
will be investigated in the next sections. 

The same effective temperature (3 of a quantum gas, in a field of randomly located oscillator wells, 
representing the Coulomb attraction between particles of the gas and localized metallic ions, was found 
in Ref. [53] by averaging the canonical density matrix of this system over the positions of the wells. A 
simple solution of the Kondo problem, in terms of T, was given in Ref. [25] . The well known puzzling 
dependence of T c on dopant concentration in high-T c superconductors, was explained also in terms of T 
in Ref. [32], as an effect arising from these interactions. 

3 Free Bose gas with the effective temperature scale T 

The fugacity of the free Bose gas (FBG) equals [55] 

z=i 1 > T " Tc ' (25) 

\ the root of g 3/2 (z) = dA 3 , T>T C , 



Tc = ( ~ 7TT ) 7— (26) 



where the BEC transition temperature T c equals 

^3/2 (l-V k 

and 

r fe / \ 1/2 M 

ff.W:=Ep \={2^(3/m) , d=— . (27) 
fe=i ' ' 

Since T(T) is a monotonously increasing function, the fugacity of the free Bose gas, with effective tem- 
perature scale T (FBET), equals 

f 1, T < f c , 

5=1 - " - (28) 

[ root of g 3/2 (z) = dA 3 , T > T c , 



with T c defined by the equality 



and 



/ \ i/ 2 

A = (27rft 2 /?/m . (30) 



The expressions for the pressure, free energy and energy of FBET are analogous to those for the FBG 

95/2(z) 

fix 



P = ^4^, (31) 
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p n,?*M + ^ s \ (32) 

P dX 

2/3d\ 3 V ; 

For the heat capacity Co and entropy So of the BGET one obtains 

15jVg 5 /2(z) i 3Ng 3/2 (z) 9i ~ 

( 4fc B /3 2 A 3 cosh 2 ( 7 /3 )T 2 ti + 2,3A 3 zd 9T ' 1 - J c, 

Co = I (34) 



15^95/2(1) 



4fc B /3 2 A 3 cosh 2 (7/3)T 2 rf' 



T<T C , 



where 



and 



§ = ^ 3 ^(/3 5l/2 (z))", (35) 

So = = iY(^co S h 2 ( 7 /3)fcBT 2 ) _1 {^|g^ +lnz}. (36) 

The FBET, with m equal to the mass of 4 He atom: m = muc = 6.64765 lCT 24 g and density of 4 He 
at T\ under saturated vapour pressure (SVP): d = dn c = 0.02197A [27], has the advantage that it 
allows to lower the BEC transition temperature below the FBG value To. Using Eq. (35), one obtains 
for m = muc, d = dn e , 

C = 4.93184008 10 8 ||{ - ^g 6/2 (g)p-V*+ 

9 , I \ / ^ , \ Q I Q ( %0 f ~\ \ "^"1 Cell 



i 53/2(5).93/2(l)(fcBro) 3/2 (/3 2 5l/2 (5)) _ }^ (37) 

if 7 is expressed in eV and &b in eV/K. The graphs of Co for various values of 7 are plotted Fig. 1. 
Experimental data Cn e (T) on heat capacity of 4 He at SVP, from Refs. [37] j [H], are also presented. The 
plot for 7 = represents heat capacity of FBG with m = tohc, d = dn e . 

Disagreement between the Co(T) plots and Cn e (T) at T < T\ is considerable. However, the Co(T) plot 
for 7 = 1.4928 10" 4 eV, with f c = 2.8 K, provides a slight improvement of the C (T) plot for 7 = (FBG). 
A further improvement, below IK, is visible in the Co{T) plot with T c — T\, 7 = 70 = 2.2276 10~ 4 cV, 
although at T\ there is no peak. 

The fraction of atoms in the momentum condensate (FMC) and zero latent heat, accompanying the 
A-transition in 4 He, are also more conveniently described in terms of the FBET than FBG. The FMC for 
the FBG equals 

/ T \ 3/2 

&(T) = !-(-) . (38) 



The corresponding quantity for the BGET therefore takes the form 



/ B \ 3 / 2 

^ = 1 -(||)) ' (39) 
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T[K] 



Figure 1: Heat capacity Co of FBET with m = mn e , d = dn c = 0.02197A , under varying 7. 7 = 
corresponds to FBG with absolute temperature scale T. The points are experimental results from Rcfs. 

H3, m- 
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where /3 C = P(j,f c ). The values of £ (T = 0) are given in Table 1. For 7 > 0, FMC satisfies | (0) < 1, 
similarly as in 4 He, where FMC is estimated to be equal 0.07 — 0.09 [TO], [IT] . 

There is no latent heat accompanying the A-transition in 4 He, whereas BEC of the FBG is a 1-st 
order transition and the accompanying latent heat per atom equals [26| 

T 5 35/2(1), _ ... 
L = - . k B T . (40) 

2 93/2(1) 



Huang derives Eq. (40), using the Clapeyron equation 

dP L 



dT TAv 

and the fact the specific volume of the BE condensate of FBG equals zero, implying 



(41) 



33/2(1) 

Analogously, for the FBET, 



A 3 

Av = v c = — . (42) 



33/2(1) 

and the vapour pressure at T < T c equals 

55/2(1) 



A 3 

Av = v c = — (43) 



P = ^44^. (44) 
/3A 3 



(45) 



Hence, 

dPp = 5 35/2(1) _ 

dT ~ 2fc B T2/32cosh 2 ( 7 ^i} c 33 /2 (l) ~ Tv c ' 
The latent heat per atom, at the transition of the BGET, therefore equals 

g5/ 2 2(1) ■ . (46) 

2 fc B T c /3(T c )2 cosh 2 ( 7 /(fc B T c ))33/2(l) 

Furthermore, Lq is correctly related to the entropy difference ASq between the gas phase and the momen- 
tum condensate. The entropy of the latter equals zero, hence ASo{T) = Sq(T). Thus, L = T c ASq/N. 
The decreasing values of Lq, for increasing 7, are given in Table 1. 
Table 1 

f c [K] 7 [eV] £ (0) L [eV] 

3.1486 1 L = 3.448585 10- 4 (47) 

2.8 1.4928 10" 4 0.59186 2.7306 10" 4 

2.17 2.2276 10~ 4 0.23187 1.4982 10" 4 
The above discussion shows that the FBET has several advantages over the FBG, as a model of low 
temperature 4 He thermodynamics. However, adjustment of T c to T\ requires setting 7 = 70, which leads 
to a flattened Co (T) graph without any peak. Further modification of the FBET is therefore necessary. 



9 



4 Bogoliubov-Huang approach to hard-sphere Bose gas with T 
scale 

The substitution Vn e —> Vi(ro = 0), performed in Section 2, amounts to considering only the attractive 
part of Vhc hi a close vicinity of the minimum. Furthermore, the thermal properties of a gas with such 
attraction are the same as those of a free gas with the temperature scale T, the FBET, if measurements 
are restricted to 1-particle obscrvablcs. Discrepancies between these properties and those of 4 He show 
that the 2-body interaction considered so far is oversimplified and, therefore, that the repulsive character 
of Vhc, at small distances, should be also accounted for by theory. These observations lead naturally to 
the concept of a hard-sphere Bose gas with the effective temperature scale T (HSET). 

The thermodynamics of the hard-sphere Bose gas (HSBG) was studied by Lee et al. [5] - [8]. The 
Hamiltonian has the form 

H = J2( Sk ~ ^) a k«k + ^jX! a p a q%>+ka q -k, 

k ' ' k p,q 

where e k = h 2 k 2 /2m, g = Airah 2 /m, a denoting the hard-sphere diameter and 

£/(*):= £[/(*)- 4 (48) 

k k 

if f(k) = Ak~ 2 + 0(k~( 3+r >>), rj > 0, as k — > oo. According to Bogoliubov [3j, [4], the interaction terms 
in H, crucial for the description of BEC, are those containing a 2 ,, a^ 2 or a^ao, owing to the observed 
macroscopic occupation of the zero-momentum 1-particlc state. By cancelling the remaining terms, 
Bogoliubov obtains the following reduced Hamiltonian : 

H r = ^(e/,- - M)a£a k + r|| a o a o Y a k a k 

k ' ' k#0 

+ 2|A| 5Z[ a o a o( a k a k + a -k a -k) + a k a -k a o + a S 2a -kflk] + ^| a o 2a o- ( 49 ) 



Owing to the macroscopic occupation of the k = mode at low temperatures, the operators ao/y |A| 
and ag/-\/|A| in H T can be replaced by complex numbers c, c, respectively. In this manner, Bogoliubov 
arrives at the Hamiltonian (see Ref. [9]) 

H B = Yj( £k _ ^ + .°l C | 2 ) a k a k + 
k^O 

+ \dYl H c l 2a k a k + |c|V k a_ k + c 2 a k al k + c 2 a_ k a k ] + i<?|c| 4 |A| - m|c| 2 |A|. (50) 

k#0 

The gauge transformation a k — > a k e* v , with adjusted ip, maps c— > |c|. Let us denote |c| 2 — > £d, where 
according to Bogoliubov [3], £ is the fraction of all atoms belonging to the momentum condensate. 
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The canonical transformation 



a k = u\J>k + Vkb*L k , u^-v^ = l, (51) 

diagonalizes Hb to the form 

Hib = Yl E ^K^ + ^ k - A*) - ^ d \ K \ + \dt?dN, (52) 

where 

hk = e k -n + 2£gd, (53) 



Eik = ^fi k -hl, h k =tgd, (54) 

and 

»H(£; + 1 )< "H^" 1 )- (55) 

The theory of BEC (e.g. Refs. [17] . [29] ) shows that the ground state of a free Bose gas becomes 
macroscopically occupied when 

lim/iA = eo as |A| —> oo, (56) 
(fi\ denoting the solution of the equation 

|A|- 1 Trzcxp(-/3i/ )(l - z cxp(-/3i/ )) _1 = d 

and the eo the lowest eigenvalue of Ho). The corresponding condition for fj,, in the BEC phase of -Hib, is 
thus 

inf lim E lk (fi) = lim E lk (fj, Q ) = 0, (57) 
m fc— >o fc— >o 

under the restriction, 

^ife(Mo) > 0, for all k g R 3 . (58) 

Eq. (56) reduces to 

fa-&d)(/i-3&d) = 0, (59) 

and only the solution //q = ^<?rf satisfies E'ifc > 0, for all k. In the BEC phase, the system may be 
therefore cquivalently described in terms of 



H B = E kKbk + ( Ek - M - l^gdN, (60) 

with (j, = and 



2' 

k#0 k^O 



f k (^d,a) = e k +(,gd 1 E k (£,d,a) = yj f 2 k - h\, (61) 
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whereas in the gaseous pase, Hb = Hq (Refs. [3], [7], [9]). 
According to the definition (48), 



^ ^[fcv^T^!- fc 2 - *2 + (63) 



where fc 2 = 8irat;d. Passing to the limit |A| —> oo, one obtains [5] 

1^' . . h 2 \K\kl f°° , 2 , ^ 2 1 , y2S 2 |A|/c5 

£ fc - A- = / <% 2 ?/vV + 2 - y 2 - l + —t = V 1C . ' 2 ' a , 64 



which yields 

H B = NE (Z, d, a) + |- V fcv/fc 2 + 2fc2fe*6 k , (65) 



2m 

k^O 



with 



£b«, d, a) = 2nad^[^J^ - (66) 
m 15 V 7r 

Lee and Yang in Ref. [3] emphasize the correct linear dependence on k of the excitation spectrum 
of Hb, as k — > 0, consistent with the presence of low-energy phonon excitations in liquid 4 He. As for 
they do not apply Eq. (56) in the range T < T c in Ref. [S] and, as a consequence, the scalar term of the 
transformed HSBG Hamiltonian is different in their theory. 

The approach to HSBG thermodynamics, developed in Refs. [S] - [8], which consists in minimizing 
the free energy 

FinfoP) = -/T^n Trexp[-Mi B ], 

(constrained by the equation for z = exp(/3/i)) with respect to £, provides only qualitative agreement of 
the resulting heat capacity with experimental data on He. 

In the present formulation, the free energy of N atoms in the BEC phase equals 

an r°° 

F B (P,Z,d,a) = NEo(Z,d,a) + dpp 2 ln(l - exp[-£ p (£, d, a)]), (67) 

y/nX^pd Jo 

where 

- 2 = 2ttM ft ^ d fl) = Jpi+^X^ad. (68) 

TO v 

In this preliminary investigation, we adopt the method presented in Ref. [7] (Chapter 3), with the 
following approximation .Fb(/9, £o, d, a) to the free energy: 

F B 0,£ o (P),d,a) := F B 0,£ o 0),d,a). (69) 

A. The energy and heat capacity 
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T[K] 



Figure 2: The heat capacity of HSET (Eq. (72)) with m = to Ho , d = d Kc = 0.02197 A~ 3 (i.e., under 
SVP at T\), 7 = 2.2276 10~ 4 eV, a = 3 A. The points are experimental results, for varying 4 He density 
at SVP, from Ref. [27j. 
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The energy U of one gram of He (expressed in calories) then equals, 

U(fi, io, d, a) = 5.757230496 1Q 3 ^ ^B^,|o,d,a) (7Q) 

Op 

(for j3 and Fb expressed in eV), where 

NE (£o,d, a) - f°° dpp 2 ln(l - cxp(-£ p )) 



9/3 v^\ 3 pdJ 

8Na§o r°° dp^ 



p 2 +4A 2 < d(exp(£; p ) - 1) 



The heat capacity of this sample obtains by numerical differentiation of U w.r.t. T: 

CyiT) = €MM*>. (72) 

Obviously, in _Fb and C/b we substitute d = dn e and 

£ (^) = 1-(S^) 3/2 , for T <T\. (73) 

For T > T c = T\, £ = 0, hence, = -Ho 13- I n the gaseous phase the energy is thus given by Eq. (33). 
By restricting the range of a to values consistent with the graph of the interaction potential Vn e [15] . 
we find the resulting CyiT) plot to be best fitting for a = ciq = 3 A to experimental data on 4 He heat 
capacity measured at SVP. It is shown in Fig. 2. Cy(T) exhibits a peak at T\, which is much smaller than 
the one observed experimentally, but below 1.8 K, Cy(T) agrees with experiment up to a small error. A 
comparison of Cy(T) with 4 He heat capacity, measured at constant density d^ c = 0.0225825 A 3 (Ref. 
[27] ) reveals similar deviation from experiment. Another point worth emphasizing is the appearance of 
singular integrals in the analytic expression for Cy(T), which result by differentiating the last term on the 
r.h.s. of Eq. (71) w.r.t. T . In the limit T — > T\, the resulting expression contains a linear combination 
(with coefficients linear in £o ) of two divergent integrals, viz., lim/(a;) and lim A(x), as x — > 0, where 

dpp 5 



Ix) 



K{x) 



o (exp(p^p 2 + x) - l)(p 2 + a;) 3 / 2 ' 



dpp A cxp(p-\/p 2 + x) 
o (exp(p^p 2 + x) - l) 2 (p 2 + x) 



We find that 



Thus, 



Inn = Inn = 1. (74) 

x— S-0 7T x— S-0 7T 



hm eo(/8)(/(eo(/3))+A'(e (/3))) - 0, (75) 

T-»T.~ 
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Figure 3: The plots of the momentum distribution function n(p, £,o(T)) for the FBET and HSET with 
a = 3AatT=lK. The points arc experimental results from Rcf. [30 . 

so CV(T) does not diverge as T — >• . It is well known that the heat capacity of 4 He, measured at 
constant volume, is finite at T\ [28j . but the experimental plot is best modelled by a function singular at 
T\ [27], [28]. The Cy(T) plot, which obtains by applying the Angelescu-Verbuere-Zagrebnov refinement 
of Bogoliubov's approximation [5], is almost exactly the same as the one in Fig. 2. 

B. Momentum distribution, momentum condensate and normal fluid density 
Theoretical plots of the momentum distribution (MD), fraction of momentum condensate (FMC) and 
normal fluid density (NFD) agree with experiment only qualitatively. The plots of normalized MD: 

Tr 



n(p,£oCO) = - t ——2 r , (76) 



[ d 3 kTr( a k ak 

for T = 1 K, 2.27 K arc shown in Figs. 3 and 4. The FMC, equal 

£(T, fo(T)) = 1 - 1 £ Tr r^^,..., : = 
iV ^ cxp[/3£ fc (£ , d, a)^ 6 k] - 1 

= 1--L^ r dpp *\ fJklM coth^%^ - l}, (77) 
^/ndX 3 Jo tE p (£ ,d,a) 2 J 

where f p (£o,d,a) — p 2 + 2a| A 2 d, is plotted in Fig. 5 for 7 = 70, a = ao, together with £o(T) for the 
same 7. Fig. 6 shows that inclusion of the hard-sphere interaction and treatment of this interaction 
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Figure 4: The plot of the momentum distribution function n(p, 0) for the FBET and HSET with a = 3 A 
at T = 2.27 K. The points are experimental results from Ref. [3"U] . 




Figure 5: The FMC £(T,f (T)) of FBET (a = 0) and HSET for a = 3 A. For 4 Hc at T = OK this 
fraction is estimated to be between 0.07 and 0.09 [TO] . 
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Figure 6: The fraction of HSET normal fluid (Eq. (78)) for d = g?hc, m = rnnei a = 3 A. In 4 He at 
T = K this fraction equals zero. 

by the Bogoliubov method, spoils the FMC graph: £(0) = 0.5224, whereas £ (0) = 0.2319. However, 
the overall description of Hell is improved (compared to FBET) if, apart from heat capacity, one also 
takes into account d n (T)/d, where d n (T) is the normal fluid density (NDF). In Hell, d n (0) = and 
£(0) + d n (0)/d ~ 0.09. d n for an interacting Bose gas, in the Bogoliubov approximation, was computed 
in Ref. [ST]. The result is analogous to NFD for the FBG: 

8 f°° - - - - 

d n {T,i (T)) = ^ dpp 4 {e E ^oAa) _ + (^(fo^o) _ i)-2}_ (78) 



30FA 3 Jo 

dn(T,io(T))/d is plotted in Fig. 7. One finds, d n (0, £o(0))/d = 0.1896 for 7 = 7o , a = a . Thus 

e(0) + ^|l = o.712, (79) 
a 

which is an improvement, as regards 4 He thermodynamics, of the FBG result 

fo(0) + ^M = l. (80) 

The plot of d n (T, £o(T))/d, for 4 He values of d, m and a = 3 A, is depicted in Fig. 6. 
C. Latent heat of the transition 

The pressure of the HSET in the BEC phase, as implied by Eq. (67), equals 

4 f°° 

P = -^stW dqq 2 \n(l-exp[-E q })-dE (^d,a). (81) 
y/irpX 6 Jo 

The Clapeyron equation (41), for P given by Eq. (81), yields the same latent heat of the BEC transition 
for the HSET as for the FBET. 
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5 Hard-sphere Bose gas with dressed condensate bosons 



Creation and annihilation processes of atoms in the momentum condensate of a Bose gas are represented 
by the operators Oq, ao, which appear, in the Bogoliubov-Huang theory, on the same footing as those 
describing atoms with nonzero momenta. However, owing to their immobility and the small fraction they 
constitute (e.g. in helium II), the condensate atoms can be expected to be particularly predisposed to 
attract the excited ones and thus act as centres surrounded by uncondensed atoms. From such point of 
view, it would be more appropriate to treat the condensate bosons as atoms surrounded by excited bosons. 
This can be achieved by performing a unitary transformation ao — > ao in H, with ao, aj representing 
the dressed bosons in the condensate. This type of procedure was applied by Frohlich in Ref. [33] to the 
electron-phonon Hamiltonian H c \_ p h. The physically meaningful terms of the transformed -ff c i_ p h were 
incorporated by Bardeen, Cooper and Schricffcr into their theory of superconductivity [34] , 
We consider the simple displacement transformation 

a = a + u, (82) 

with ^eR 1 representing the cloud of excited bosons accompanying the immobile one, as a possible trial 
realization of the mapping ao — > ao. Glauber demonstrated in Ref. |35j that the mapping (82) can be 
accomplished by the unitary transformation 

D{v)- 1 < M) D{v)=a + v, (83) 

where 

D(u) = exp[u(a* - a )]. (84) 

In the transformed Hamiltonian Hid = D{y)~ 1 H r D(y), we perform Bogoliubov's substitution ao — > 
\/£d|A| and, furthermore, v — s- y/i]d\A\/2, where rj/4 > is the fraction of bosons which participate in 
the dressing of those in the momentum condensate. Since < £ <C 1 in 4 He, therefore 77/4 can be also 
expected to be small. Thus ^/f > £, y/v/^ > v/^> so the terms of the reduced interaction in Hid, most 
significant for the description of BEC in 4 Hc, arc those linear in yf^q and £. The latter were considered 
in Sec. 4. In this section, we take into account only the former ones and disregard all others. In this 
manner, we obtain the Hamiltonian 

^9\/Cnd^2 (ak a k + al k a-k + ak a -k + a -kOk) - Mv^? rf l A l- ( 85 ) 

k#0 

It should be noted at this point that Bogoliubov rejected in his paper [3] all terms of the 2-body 
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interaction linear in ao, cLq, such as 

^Xj a o a q a p a q-q' Wm ^ a P -q a q a p a o- ( 86 ) 

p-q/0 p-q^O 

These terms cannot be included into Bogolibov's diagonalization procedure. Processes described by the 
expressions (86), representing the creation or annihilation of a single boson in the momentum condensate, 
as a result of 2-body interactions, should be expected to occur especially at temperatures T satisfying 
< 1 — T/T\ <C 1, where the number of such atoms is macroscopic, but still very small. Hjj can be 
viewed as a Hamiltonian which simulates such processes. 

The transformation (51) diagonalizes Ho and the same rearrangements , as applied in Sec. 4 to Hb, 
transform Hjj to the form 

H D = NE D & d, a D ) + |- X>\A> 2 + 2k D 6 p 6 P' ( 8? ) 



2m 

p^O 



where 



H 2 /128 a 3 £ 5 / 2 d \ . r- 

E D (£, d, a D ) = 2ira D d — (-^\/ — 2£,Vv ' k D = 8n V £ and, a D = ay/rj. 

m \ 15 V 7r / 

In the sequel we put 77 = 1, in order to restore the unsealed hard-sphere diameter a. The system 
described by H D , with 77 = 1, will be called the D-transformed HSET (DHSET). 
A. The energy and heat capacity 

For ?7=1, the energy of 1 gram of 4 He (expressed in calories) under the same approximation £ = £0, 
as in Sec. 4, equals 

U(0, f 1/2 , d, a) = 5.757230496 10 3 (^ (^ /2 , d, a) - / dpp 2 ln(l - cxp(-E p (^ /2 , d, a))) 

^-KA^pd Jo 

Wa^ 2 r dPP" ^ (89) 



^ Jo yft + AaX 2 ^ 2 d (exp(^(| 1/2 , d, a)) - 1) 

(for /3, A 2 expressed in eV). The square root t^oifi) 1 ^ 2 appearing in the coefficient of the second integral 
on the r.h.s. leads to the singularity 

of heat capacity, as T — > Tj~. The best fitting to experiment heat capacity graph, obtained for d = 
dft e ,m = rriHe and 7 = 70, a = 2.45A, , by differentiating U(j3, io(P) 1 ^ 2 , d, a) numerically w.r.t. T, 

uCPMP) 1,2 ,d,a) 

LV{1) = -jrp ) ( yi ) 

is depicted in Fig. 7. There is good quantitative agreement with experimental data on 4 He heat capacity 
at SVP below 2.1 K. For T £ [2.1 K, 2.17K] the C v , given by Eq. (91), exceeds the experimental values. 
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T[K] 



Figure 7: The heat capacity CV(T), given by Eq. (91), for d = dn e , m = ?tihc, a = 2.45 A, 7 = 70. The 
points are experimental heat capacity results from Rcf. [2 7) . 

It is therefore indeed the temperature range < l — T/T\ <C 1, where the CV(T) plot in Fig. 2 is modified 
by the Hamiltonian Hjj most significantly. Furthermore, the best fitting value a = 2.45 A for Eq. (91), 
agrees with smaller error, than 00 in Sec. 4, with the empirically established hard-sphere diameter a f=s 
2.6 A of 4 He atoms [H]. 

The latent heat of the Hd BEC transition is the same as for Hb- 
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Figure 8: The plots of the MD function n(p, £o(l) 1/2 ) for the FBET and DHSET with a = 2.45 A, 7 = 7o 
at T = 1 K. The points are experimental results from Ref. [30] 

B. Momentum distribution, momentum condensate and normal fluid density 

The DHSET Hamiltonian improves the theoretical plots of MD, PMC and NFD obtained in Sec. 4. 
These quantities for the DHSET, equal n(p, |o(T) 1/2 ), £(T, fo(r) 1/2 ), d«(T, £ C r ) 1 < /2 ), respectively, arc 
plotted in Figs. 8, 9, 10. For T = 2.27 K, the plot of n(p, | (2.27K) 1 / 2 ) is the same as in Fig. 4. 

These results, as well as the DHSET heat capacity plot in Fig. 7, show that the DHSET descrip- 
tion of helium II thermodynamics considerably improves the results of HSET theory. In particular, 
£(0,|o(0) 1/2 ) = 0.4266, d„(0,Co(0) 1/2 )/d = 0.1178 and 

COUo(O) 1 / 2 ) + d " ( °'y )1/2) = 0.5444, (92) 
a 

which is an improvement of the corresponding HSET result (Eq. (79)) by 23.5%. 
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Figure 9: The FMC function £(T, ^(T) 1 ^ 2 ) of DHSET for a = 2.45 A 7 = 7o , d = d Rc , m = m Hc - For 
4 He at T = K this fraction is estimated to be between 0.07 and 0.09. [TU] 
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Figure 10: The fraction of DHSET normal fluid density d n (T, ^{Tf/^/d for a = 2.451, 7 = 70, d = 
due, "m = mjie- In 4 He at T = K this fraction is equal zero. 
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6 Concluding remarks 

The DHSET has revealed good quantitative agreement of heat capacity with experimental 4 He mea- 
surements below 2.1 K and considerable improvement of the HSET heat capacity graph, although the 
sharp, but continuous, decrease of heat capacity observed above T\ in 4 He, is not reproduced the HSET 
nor DHSET. Other thermal properties of the DHSET also improve the HSET results, but agree with 
experimental data on 4 He only qualitatively. Further improvement of this theory is therefore necessary. 
One possible extension could include a more general form of effective temperature, another inclusion of 
all terms of H\r>, linear both in s/£rj and £, into the theory. 
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